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We investigate gauge dependence of a nonrelativistic boundstate formalism used 
in contemporary calculations. It is known that the effective Hamiltonian of the 
boundstate system depends on the choice of gauge. We obtain the gauge trans- 
formation charge of the Hamiltonian, by which gauge independence of the mass 
spectrum and gauge dependences of the boundstate wave functions are dictated. 
We raise two questions of practical and physical interest, and provide answers to 
them. 

1 Introduction 

Recently there has been much progress in our theoretical understanding of 
nonrelativistic QED and QCD (NRQED and NRQCDL hqundstatcs such as 
positronium, T and remnant of toponium boundstatcs EJlij. A notable char- 
acteristic in these new developments is that the conventional Bethe-Salpctcr 
equation is no longer being used to calculate the spectrum and wave func- 
tions of boundstates. Instead, one starts from the non-relativistic Schrodinger 
equation (of quantum mechanics) with the Coulomb potential. Then one adds 
to the nonrelativistic Hamiltonian relativistic corrections and radiative correc- 
tions as perturbations to obtain an effective Hamiltonian (quantum mechanical 
operator) valid up to a necessary order of perturbative expansion. Effective 
Hamiltonians used in these new formalisms are known to be dependent on the 
choice of gauge. 

We report here our recent achievementsu-H on investigations of gauge de- 
pendence of a nonrelativistic boundstate formalism used in contemporary cal- 
culations. Our motivations for the study are: (I) In the present frontier calcu- 
lations of higher order corrections to physical quantities of boundstates, often 
the Feynman gauge is used to calculate typically ultraviolet radiative correc- 
tions whereas the Coulomb gauge is used to calculate corrections originating 
typically from infrared regions. Therefore, it is desirable to clarify gauge de- 
pendences of the formalisms actually used in these calculations. (II) We would 
like to find transformations of boundstate wave functions when we change the 
gauge-fixing condition. We may apply these transformations to study vari- 

a Talk given at the International Symposium on "Quantum Chromodynamics and Color 
Confinement" (Confinement 2000), Osaka, Japan, March, 2000. 



1 



ous amplitudes involving boundstates. Since a physical amplitude is gauge 
independent, once we know how the wave function transforms, we know how 
other parts of the amplitude should transform to cancel gauge dependence as 
a whole. This would provide a useful cross check for identifying all the con- 
tributions that have to be taken into account at a given order of perturbative 
expansion. 

A formal argument of gauge dependence/independence of boundstate for- 
malisms goes as follows. It is well known that a quark-antiquark boundstate 
contributes a pole to the quark-antiquark 4-point Green function: 



SI/ 7->\ ^ ^IA.P 

G{pqP) = — — ' — •— + (regular as P° - u p), (1) 



An infinitesimal deformation of the gauge-fixing function in the QCD La- 
grangian, F — > F + SF, induces a variation 



Any physical amplitude (/; out|i; in) which involves the quark-antiquark bound- 
state contributions includes the above Green function G(pqP) as a part of it. 
Since the initial and final states satisfy the physical state conditions Qb | i; in ) = 
Qb I /;in) = and the theory is BRST invariant, the amplitude is gauge in- 
dependent. Hence, the boundstate poles included in the amplitude are also 
gauge independent. 

Nevertheless we may raise two intriguing questions: (I) Does the Green 
function G(pqP) include any unphysical pole, which does not contribute to the 
physical amplitude, close to or degenerate with one of the physical boundstate 
poles? A typical example is the i?^-gauge for electroweak interaction where an 
unphysical pole (k 2 — £,M^ +ie)~ 1 is included in the gauge boson propagator. 
(II) Is the boundstate wave function a physical observable? The answer is 
obviously no since the wave function is gauge dependent as a consequence of 
the gauge dependence of the effective Hamiltonian. Then, what is the (gauge- 
independent) physical quantity which can be thought as a counterpart of a 
boundstate wave function? 

We will answer to these questions below. 

2 Effective Hamiltonian and Transformation Charge 

Let us sketch the outline of our argument. The quantum mechanical Hamilto- 
nian is determined from perturbative QCD order by order in expansion in 1/c 
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(inverse of the speed of light): 



H = H Q + - H x + \ H 2 + ■ ■ ■ . (3) 



c cr 

Since quark and antiquark inside a heavy quarkonium are non-relativistic, the 
expansion in 1/c leads to a reasonable systematic approximation. Presently the 
Hamiltonian is known up to C(l/c 2 ). There exist several different definitions 
of an effective Hamiltonian for the NRQCD boundstates beyond leading order. 
We introduce an effective Hamiltonian defined naturally in the context of time- 
ordered (or "old-fashioned") perturbation theory of QCD. Then we obtain a 
transformation charge Q (quantum mechanical operator) such that the effective 
Hamiltonian and the boundstate wave function change as 

5H cS (P°) = [H eS (P°) - P°] iQ(P°) - iQ\P°) [H eS (P°) - P°] , 
Sip = —i Q ■ ip 

when the gauge-fixing condition is varied infinitesimally. Written in terms of 
the BRST charge and the field operators in the QCD Lagrangian, Q is given 
by 

(p,-p,\,\\Q(P°)\q,-q,\\\') 

A", A" 

xg-i(q',q;\",\",\',\';P°). 

Gauge independence of the spectrum can be shown using the transformation. 
Since Q has no pole ~ (E — M v + ie)^ 1 , there is no unphysical state which 
contributes a pole to the Green function Q that is degenerate with or close 
to one of the poles of the physical boundstates of our interest. Stating more 
explicitly, there is no unphysical boundstate with a binding energy ~ a 2 s m. 

3 Application 

It is known that the top quark momentum distribution in the ti threshold 
region at leading order is proportional to the absolute square of the wave 
functions of (would-be) toponium boundstates in momentum space E3. The 
momentum distribution will be measured at future collider experiments, so 
we wilL-be able to probe the boundstate wave functions. As we have shown 
in Reftil, wave functions of boundstates are gauge dependent beyond lead- 
ing order. We verified that this gauge dependence is cancelled by that of the 
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final-state interaction diagrams at C(l/c) when calculating the top quark mo- 
mentum distribution. In other words, a boundstate wave function mixes with 
the final-state interaction diagrams by gauge transformation. 

Physically b quark emitted in top decays carries the color charge and in- 
evitably interacts with the gluons which were also responsible for the bound- 
state formation. Since the top quark momentum is reconstructed from the 
momenta of its offsprings, it is natural that the QCD interaction of the off- 
springs mixes via gauge transformation in the determination of top momentum 
distribution. ._. 

At the same time, this shows that the present calculations^! of the top 
momentum distribution at C(l/c 2 ) are gauge dependent, i.e. they vary if we 
transform the gauge infinitesimally from the Coulomb gauge, since they do 
not include the final-state interaction diagrams. The example at 0(1/ c) sug- 
gests how gauge cancellations should take place in the complete amplitude at 
C(l/c 2 ) which has not been obtained yet. 

4 Summary 

• We used the BRST symmetry to formulate our arguments, which enabled 
us to discuss gauge dependence of the NRQCD boundstate formalism 
rigorously. 

• Presently, there exist several different definitions of an effective Hamil- 
tonian beyond leading order. We introduced an effective Hamiltonian 
defined naturally in the context of time-ordered (or "old-fashioned") 
perturbation theory of QED/QCD. Then we obtained a transformation 
charge Q (quantum mechanical operator) when the gauge-fixing condi- 
tion is varied infinitesimally. Also, gauge independence of the spectrum 
is shown using the transformation. 

• For illustration: (l)we calculated the transformation charge Q at next-to- 
leading order; (2)we demonstrated gauge cancellations among diagrams 
by examining an infinitesimal gauge transformation of the amplitude for 
a qq boundstate decaying into q'q"W + W~ . From the latter example, 
one can deduce that the present calculations of the top momentum dis- 
tribution in the tt threshold region at next-to-next-to-leading order are 
gauge dependent. 
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